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Abstract 

Single molecule DNA experiments often generate data from force versus extension 
measurements involving the tethering of a microsphere to one end of a single DNA 
molecule while the other is attached to a substrate. We show that the persistence 
length of single DNA molecules can also be measured based on the recoil dynamics 
of these DNA-microsphere complexes if appropriate corrections are made to the 
friction coefficient of the microsphere in the vicinity of the substrate. Compari- 
son between computer simulated recoil curves, generated from the corresponding 
Langevin equation, and experimental recoils is used to assure the validity of data 
analysis. 
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1 Introduction 



One of the most important characteristics of the DNA molecule is its high 
mechanical flexibility. Thanks to this, the entire genome of each living organ- 
ism, which can reach several meters in length, fits inside each of its cells. The 
flexibility of the DNA molecule also plays a key role in all its cellular func- 
tions, such as replication, transcription and recombination. With advances in 
micromanipulation techniques during the last decade single molecule DNA 
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experiments have become common, and are explored through many different 
approaches [1]. 

Once free in solution, single DNA molecules present entropic elastic behavior, 
well described, in both low and high stretching regimes, by the worm-like chain 
(WLC) model [2-4]. The flexibility of the polymeric chain can be inferred 
through the ratio A/L, where A is its persistence length and L its contour 
length. Smaller A/L ratios correspond to more flexible chains. For the DNA 
molecule, a complete chain can have contour length L varying from around 
2fxm to more than lm, while the persistence length A is of the order of 40 — 
50nm or 120 — 150 base pairs. 

The usual method to obtain A employing an optical trap [5] is based on force 
versus extension measurements [6-9]. One end of a DNA molecule is attached 
to a coverglass surface while the other end is tethered to a polystyrene micro- 
sphere (Fig. 1). The optical trap can then be used to capture and move the 
microsphere, stretching the DNA molecule. Obtaining force versus extension 
data for a DNA molecule, however, requires a very precise characterization 
and calibration of the trap, allowing the determination of the trapping force 
as a function of the position of the trapped microsphere. Here we examine 
an alternate method to determine the persistence length A: once the DNA is 
stretched, one removes the trap and follows the motion of the microsphere in 
time, using videomicroscopy [8]. 

Using this method, Feingold [10] obtained a value of 152nm for the persistence 
length of DNA, approximately three times larger than the value obtained 
from force versus extension measurements. The proposed explanation for this 
discrepancy was that during the initial part of the recoil, when the DNA is 
more than 80% stretched, the motion cannot be considered quasistatic, such 
that the force at each position is time dependent. Recently, Bohbot-Raviv 
et al [11] considered a complete nonequilibrium theory for the relaxation of 
highly stretched semiflexible polymers like DNA and found a better agreement 
with Feingold's experiment. 

We have performed recoil experiments in which the friction between the micro- 
sphere and the coverglass is carefully taken into account and found persistence 
lengths comparable to the values obtained from force versus extension mea- 
surements just using the standard quasistatic approximation, without the need 
of a nonequilibrium theory for the relaxation of the DNA. We also perform 
computer simulations of the recoil using the corresponding Langevin equation 
to assure the validity of the assumptions used in the data analysis. 

The balance of this paper is organized as follows. In section 2 we describe 
the recoil method and the motion of the microsphere. Section 3 details the 
simulation of the recoil dynamics of the microsphere. The procedures adopted 
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Fig. 1. Microsphere anchored to a coverglass surface by a single tethered DNA 
molecule. 

to obtain experimental recoil curves are described in section 4. In section 5 
we describe the analysis of simulated and experimental recoils leading to our 
results. Conclusions and final remarks are presented in section 6. 



2 Experimental Technique 



In essence the recoil technique introduced by Shivashankar et al [8] is very 
simple. As in the force-based experiments, optical tweezers are used to find 
and trap a microsphere anchored to a coverglass surface through a single teth- 
ered DNA molecule. The trapped microsphere is moved, stretching the DNA 
molecule. The tweezers are then cut off suddenly, releasing the microsphere, 
which is dragged through the fluid by the recoiling DNA molecule. Following 
the recoil of the microsphere using videomicroscopy, we obtain a recoil curve 
R(t), from which we can extract information about A and L. 

According to [8], the recoil movement can be considered quasistatic since the 
equilibration time for the DNA molecule is very short compared with the time 
scale for the motion of the microsphere. This means that, during relaxation, 
the DNA essentially passes through a sequence of equilibrium stretched states 
[12]. Because the time resolution of our videomicroscopy experiments ((l/30)s) 
is large compared to the ratio m/7 (ss 10~ 7 s) for the microsphere, inertial 
effects are negligible in the analysis of the recoils. As a result, the average 
motion of the microsphere can be described by 

7^ = -F dna {r\ (1) 



where (R) is the average position of the microsphere, r is the distance between 
the extremities of the DNA molecule, 7 is the Stokes friction coefficient and 
Fdna{f) is the entropic elastic force for the DNA molecule given by the WLC 
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model [2-4], 



F d na(r) 



k b T 



A 




) 



(2) 



where the final term represents the correction introduced by Bouchiat et al in 
[4] with a 2 = -0.5164228, a 3 = -2.737418, a 4 = 16.07497, a 5 = -38.87607, 
= 39.49944 and 0:7 = —14.17718. The bouyancy force on the microsphere 
is 5fN, negligible compared to F dna , which is of the order of pN during most 
of the recoil movement. 

Far from any boundary, the Stokes friction coefficient between the microsphere 
and the surrounding medium is 

7 = 67ir]a, (3) 

where 77 is the viscosity of the medium and a is the radius of the microsphere. 
There are, however, important considerations to be made regarding the use 
of Eq. (1). The wall-drag effect caused by the coverglass nearby changes the 
friction coefficient 7 of the microsphere [13]. As a result, 7 can no longer be 
considered constant and isotropic, and must continuously be corrected as the 
microsphere moves. This correction depends on the height h of the microsphere 
with respect to the coverglass below. 

We can integrate Eq. (1) directly. However, we cannot obtain a simple explicit 
analytic expression for R(t), forcing us to devise an alternative procedure to 
fit the experimental data. We chose to apply a discrete derivative to the ex- 
perimental recoil curves, generating data that may be fit directly with Eq. (1). 
Since the microsphere undergoes Brownian motion while recoiling, the numer- 
ical derivative will, however, be noisy. We therefore smooth the recoil curve, 
minimizing the Brownian noise, prior to performing the numerical derivative. 
To test the validity and effects of the adopted procedure, we compare our 
result with the simulated the recoil dynamics of a tethered microsphere. 



3 Simulation 

To simulate the recoil dynamics of the microsphere we solve, numerically, the 
Langevin equation for the system, using F dna as an external force [14,15]. In 
the non-inertial limit, we have 



dR 
~dt 



F dna (r) F(t) 



(4) 
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where J-{t) is the Brownian force, which has the properties {^{f)) = and 
{F{t)F{t')) = T5(t - t'), where T = 2 1 k B T. 

The discrete form of Eq. (4) is 

R(t + At ) =R(t)-^^ + ^i r ,, (5) 

7 7 



where £ n is a Gaussian random variable with (£ n ) = and (£ n £m) = Sum- 

Until this point, we have treated the motion of the microsphere using a ref- 
erence frame directed along the force axis of the DNA molecule. This one- 
dimensional reference frame is not appropriate to our purposes, since it does 
not contain the experimental reference frame which we want to simulate, in 
which the microsphere undergoes three-dimensional Brownian motion. There- 
fore, we extend Eq. (5) to three-dimensional space, obtaining 

v,, A ^ Y(A F dna (r)At X t (t) y/TAt 

Xi(t + At) = Xi(t) - 1 4„,i, (6) 

7 K 7 



with % — 1, 2, 3 and 



i=3 



— a, — R — a, 



(7) 



where the set {Xi} represents the coordinates of the center of mass of the mi- 
crosphere while r represents the relative distance between the two extremities 
of the DNA molecule. Since the DNA is attached to the surface of the micro- 
sphere, the two coordinate systems are connected through Eq. (7), where a is 
the radius of the microsphere. 



3.1 Correcting^: The wall- drag effect 



Use of a constant and isotropic Stokes friction coefficient 7 is justified if and 
only if the fluid velocity attains constant value v far from the microsphere, 
in all directions. When an obstacle interferes with the fluid flow near the 
microsphere, such as the coverglass surface in the present case, this assumption 
is no longer valid and 7 depends on the direction of motion and on the distance 
from the microsphere to the obstacle [13]. According to the analyses of Faxen 
[16], Brenner [17] and Goldman et al [18] in the early and mid 20th century, 
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for a microsphere of radius a whose center of mass is at a distance h from a 
planar surface, the friction coefficient 7 is anisotropic, with 




for motion parallel to the planar surface (xy plane), and 



7^ = Qur/a x 



4 . ~ n(n + l) 

- sinn a > — 

3 ii(2n-l)(2n + 3) 



x 



2 sinh [(2n + l)q 


+ (2n + 1) sinh 2a 


4 sinh 2 


(n + \)a 


- (2n + l) 2 sinh 2 a 



- 1 



(8b) 



for motion perpendicular to the planar surface (z axis), where a = cosh _1 (/i/a). 
The exact correction for j±, shown in Eq. (8b), can be replaced by a far simpler 
approximation, given by [17] 



7_i_ ~ Qnrja ( 1 + 



h 



ic) 



with an average error < 5% for h/a > 1. The loss of isotropy in 7 alters 
Eq. (6), which becomes 



X i(t + A t ) = X,( t )-^^^ + ^„„ ,= 1,2, (9a, 

7|| it 7|| 



^HA^A-.W-M^M,^, , 9b) 

7_l il! 7 ± 



where Eq. (9a) generates the x and y components of the recoil curve and 
Eq. (9b) the z component. This equation, however, does not take into account 
the fact that the height h of the center of the microsphere cannot be less than 
its radius a, in which case the microsphere is touching the coverglass. At this 
point, interactions between the microsphere and the coverglass can no longer 
be neglected and Eq. (8) can no longer be safely used. We have nevertheless 
observed that, after the recoil, the microsphere exhibits Brownian motion in 
the plane of the coverglass for quite some time, the same behavior being ob- 
served for most of the free microspheres as well. This observation indicates 
that the microspheres remain hovering close to the coverglass for a period of 
time before sticking to it. To determine how close to the coverglass the micro- 
spheres were hovering, we prepared a sample containing only free microspheres 
under the same experimental conditions. After letting the sample rest on the 
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Fig. 2. Average quadratic displacement as a function of time for microspheres 
hovering close to the coverglass. The movements of 20 free microspheres were 
tracked for 40s and combined into a single curve. The difusion coefficient 
D = (0.063 ±0.001)/xm 2 /s) 

was obtained by fitting this curve to a two-dimensional 

random walk. 



microscope for more than 2 hours, we measured the microsphere diffusion co- 
efficient for motion parallel to the coverglass (Fig. (2)). The measured value 
was (0.063 ± 0.001)/im 2 /s, approximately 2.75 times less than the theoretical 
free difusion coefficient, = 0.173/im 2 /s, for a = 1.42/im. Assuming this 
discrepancy is caused by the increased friction close to the coverglass and ap- 
plying Eq. (8a) we obtained an average hovering distance of approximately 
1.47 /im, indicating that the microspheres in fact were not touching the cover- 
glass, their point of nearest approach remaining approximately 0.05/im above 
it. The simulation routine was then altered to reproduce this effect, not let- 
ting h drop below its minimum value of 1.47/im during the recoil, as observed 
experimentally. 

The initial conditions of the microsphere for all simulated recoils are X\ (0) = 

J (0.85L + a) 2 — hi, X 2 (0) = and X 3 (0) = h , where we have assumed that 
the DNA molecule is at least 85% stretched, initially. The coordinate system 
is taken so as to place the direction of stretching in the xz plane. The time 
increment, At, was set to (l/30)s, matching videomicroscopy time resolution 
for better comparison between simulated and experimental curves (which also 
proved adequate for numerical convergence). A simulated recoil curve is shown 
in Fig. (3). 
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Fig. 3. A typical simulated recoil curve. The numerical solution of Eq. (9) is plotted 
for each direction. For this simulation we used ho = 5/xm, A = 45nm, L = 17/xrra 
and T = 25°C. 

4 Experiment 

4-1 Experimental Setup 

The experiments are made on a Nikon TE300 Eclipse inverted optical micro- 
scope coupled to two CCD cameras and a piezo-driven stage (Fig. 4). We use 
a SDL 5422 near infrared laser (X = 832 nm ) focused through an infinity cor- 
rected objective lens (Nikon plan apo, DIC H, WOx, 1.4 NA, oil immersion) 
to create the optical tweezers. Imaging from CCD 1 is used to locate and move 
suitable microspheres into position while imaging from CCD 2, filtered from 
the near infrared light, is videorecorded for posterior analysis. 

4-2 Sample Preparation 

Samples are prepared using the method described in [8], with slight variations. 
We prepare a solution containing PBS (150 mM Na + ) pH 6 and microspheres 
(Polybead Polystyrene 2.8/xm microspheres, p = 1. 05 g/ml,n = 1.6). The den- 
sity of microspheres in the solution is around 450/ fil. For every 395 \xl of solu- 
tion, we add 5 fil of DNA solution (New England Biolabs X DNA 553 fig /ml), 
previously heated in a thermal bath at 62°C for 5 minutes, mix gently and in- 
cubate together for 20 minutes. The solution is then pipeted into cells (~ 4mm 
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Fig. 4. Schematic of the experimental setup. (IL) light source, (LI) objective lens 
20x, (L2) objective lens lOOx, (L3) condenser, (AP) anamorphic prism, (Ml) di- 
croic mirror, (Fl) near infrared filter. 

radius, 3 mm height) on coverglasses and covered with cut pieces of coverglass 
to prevent excessive evaporation. At this pH condition, DNA molecules bind 
to the microspheres and to the coverglass, preferably by its extremities [19]. 
After 24h of incubation at room temperature, the cells are opened and gently 
washed with PBS (150 mM Na + ) pH 7.4 to remove the excess of free beads 
and reset the pH for the DNA. The cells are covered again and are then ready 
for experiments. 



4-3 DNA-Microsphere System Selection 



With the samples placed on the microscope, the optical tweezers are used 
to locate a microsphere bound to a single DNA molecule. Although many 
microspheres are attached to DNA, some judgement is needed to choose a good 
DNA-microsphere system in order to avoid microspheres attached to more 
than one molecule, molecules bound to the coverglass at multiple locations and 
other possible undesirable configurations. Verifying if the selected microsphere 
presents radial symmetry with respect to stretching and recoiling around its 
apparent resting point is the minimal required test for a good candidate. Once 
a suitable DNA-microsphere pair is selected, the initial height ho for the recoil 
is set adjusting the focal distance of the tweezers. The microsphere is trapped 
and moved to stretch the DNA. The videorecording is started and the tweezers 
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cut off, starting the recoil. 



4-4 Data Extraction 



After videorecording the recoil through CCD 2, the images are analyzed using 
the Image 1.62 [20] software package. The x and y coordinates of the center 
of mass of the microsphere are extracted from each frame and combined with 
the time index of the recording to give X±(t) and X 2 (t). 



4-5 Estimate for h(t) 



As already mentioned, the appropriate correction to the friction coefficient for 
the microsphere depends on the distance h of its center from the coverglass. 
We use the optical tweezers to set the initial height h of the microsphere. 
During the recoil, the value of h can be measured from the defocusing of the 
microsphere, a very time consuming procedure. We believe that, for the re- 
coil experiments, sufficient accuracy can be achieved by using a reasonable 
estimate for h(t) guided by computer simulations. To obtain h(t) for the ex- 
perimental recoils we used averaged simulated recoils to establish how h(t), 
which corresponds to X$(€) in the simulated recoils, was related to the horizon- 
tal component of the recoil, Xi(t). We found that, because of the anisotropy 
in 7, h{t) plotted as a function of Xi(t) does not correspond to a straight line 
trajectory pointing to the origin of the coordinate system, as would be ex- 
pected if 7|| and 7_i_ were identical. Nevertheless, we observed that the average 
h(t) can be estimated, with good agreement with the simulations (Fig. (5)), 
by 

h(t) = 0.925 + (ho - 0.925) x ^ Q) , (10) 



until h{t) reaches 1.47/im, where the experimental evidence indicates that the 
microspheres interrupt, at least for a time longer than the duration of the 
recoil, their vertical motion. During the most relevant region of the recoil, the 
distance from the bottom of the microsphere to the coverglass is larger than 
0.5/xm, a value much larger than any roughness of such a good optical quality 
coverglass used. 
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Fig. 5. The height h of the center of mass of a microsphere is plotted as a function 
of its horizontal position, (o) Average of 250 simulated recoils showing (h(t)) as a 

function of the corresponding (Xi(t)}. ( ) Estimate for h(t), given by Eq. (10). We 

have tested Eq. (10) with different values of h,Q, A, L and temperature and obtained 
similar results. Parameters: ho = 5fim, A = A5nm, L = 17/im and T = 25°C. 

5 Data Analysis 



Using the previously chosen coordinate system, in which the molecule is stretched 
along the xz plane, the recoil is primarily along the x axis. We may then use 
the x component of Eq. (1), given by 

dXi _ F dna (r) X 1 
dt 7// R ' 



to fit our data. The y component of the recoil will fluctuate around 0, allowing 
us to write 



R w J XI + h 2 . (12) 



Analysis of the recoil curves (Fig. (6)) consists of three steps. First, the re- 
coil curve is smoothed using the Stineman smoothing function from Kaleida- 
Graph [21] software package, filtering out high-frequency Brownian noise. The 
smoothed recoil curve is then used to evaluate the finite difference The 
resulting velocity data are fit according to Eq. (11). 
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Fig. 6. Simulated (a) and experimental (b) recoil curves. The parameters of the 
simulated curve were chosen close to the values obtained from the analysis of the 
experimental recoil, in order to show that the simulated recoils are a good reproduc- 
tion of the experiments. The following values were used: A = 45rem, L = 16.5/^m, 
ho = 3\xm and T = 23°C. 
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Fig. 7. Analysis of the recoils shown in Fig. (6). a) Simulated recoil, b) Experimental 
recoil. For the simulated recoils, differences between the parameters A and L used in 
the simulation and the values obtained from the analysis of the recoil are considered 
acceptable if < 5nm for A and < 0.3fim for L. 

To check the validity of the results produced by this procedure, we compared 
the values of A and L obtained from the analysis of several simulated recoils 
with the known values of A and L used as parameters in each simulation. The 
process was repeated several times, with different values of A, L, h and T 
used in the simulated recoils. In most cases, there was good agreement between 
the values of A and L used in the simulation and the values obtained from 
the analysis of the recoils generated (Fig. (7)). We have, however, observed 
that in certain cases the analysis of a simulated recoil produced results quite 
different than those obtained from the analysis of other recoils generated with 
the same parameters. We believe that such occasional discrepancies occur 
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when a slight deviation from the average recoil trajectory is enhanced when 
the recoil is smoothed and its derivative evaluated. To minimize this effect in 
the experimental analysis, we have always recorded, analysed and compared 
at least 4 recoils for each DNA-microsphere complex. 

We analysed 21 experimental recoils from 5 different DNA-microsphere sys- 
tems in 3 different samples. The average persistence length obtained from our 
data was A = (43 ± 5)nm, a value well within the range obtained from other 
techniques [6,9,22]. 



6 Conclusion 



In recent years measurements of the persistence length of single DNA molecules 
have been made through several different methods, ranging from spectroscopic 
techniques to electron microscopy. Each of these approaches involves a partic- 
ular set of assumptions, as pointed out by Wang et al. in [6], and is susceptible 
to different sources of error. 

We have shown, based on simulated and experimental data, that the recoil 
method, in which we follow the recoil dynamics of a previously stretched DNA 
molecule through the motion of a tethered microsphere, can be used to deter- 
mine the persistence length of single DNA molecules, provided that adequate 
correction of the friction for the microsphere is made. Although it is possible 
to perform the analysis of the recoil using only on the horizontal component 
of the movement, knowledge of the vertical dynamics of the microsphere is 
fundamental, since its vertical position is the critical parameter governing the 
friction correction. Both simulated and experimental data were treated us- 
ing the quasistatic response of the DNA molecule, without considering the 
far from equilibrium dynamics proposed in [11]. We observe a high degree of 
consistency between experiment and simulation. 
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